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Abstract 

The energetic exclusive two-body nonleptonic decays of B mesons are inves¬ 
tigated in the framework of the relativistic quark model within the factorization 
approximation. The heavy quark expansion is used for the calculation of form 
factors. The obtained results are in agrement with available experimental data. 
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1 INTRODUCTION 

The investigation of exclusive nonleptonic decays of B mesons represents an important and 
complicated theoretical problem. In contrast to the exclusive semileptonic decays, where 
the weak current matrix elements between meson states are involved, nonleptonic decays 
require the evaluation of hadronic matrix elements of the local four-quark operators. To 
simplify the analysis it is usually assumed that the matrix element of the current-current 
weak interaction factorizes into the product of two single current matrix elements. Thus 
the problem reduces to the calculation of the meson form factors, parametrizing the 
hadronic matrix elements of weak currents as in the case of semileptonic decays, and 
the meson decay constants, describing leptonic decays [Q. This makes the factorization 
hypothesis to be a very appealing assumption. However, strong interaction effects, such 
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as final state interactions, the rescattering of the final hadrons etc., can violate this 
approximation i,0. There are also some problems with the different renormalization 
point dependence of the initial and factorized amplitudes [fj]. Thus factorization can 
not be considered as a universal approach to nonleptonic decays. 

There were several theoretical developments which can help to justify the factoriza¬ 
tion for certain nonleptonic decays of heavy mesons. It has been shown in ref. || that 
factorization holds in the limit of large number of colours N c in QCD. The leading 1/IV C 
corrections to this limit have also been considered. Moreover intuitive arguments justi¬ 
fying factorization for the energetic nonleptonic decays were given by Bjorken || on the 
basis of the so-called colour transparency. In these decays the final hadrons are produced 
in the form of point-like colour-singlet objects with a large relative momentum. And thus 
the hadronization of the decay products occurs after they are too far away for strongly 
interacting with each other, providing the possibility to avoid final state interactions. 
Dugan and Grinstein [;3] discussed the factorization hypothesis within the heavy quark 
effective theory. They proved that factorization holds for the decays of a heavy B meson 
into a heavy D meson and a light meson, where the light quarks, which hadronize into 
a light meson, are highly energetic and collinear. Therefore we have a good theoretical 
background to expect that factorization can be applied to the consideration of energetic 
nonleptonic decays of B mesons. 

In this paper we calculate the branching ratios of the exclusive energetic nonleptonic 
decays of B mesons in the framework of the relativistic quark model on the basis of 
factorization. The heavy-to-heavy hadronic form factors, appearing in the factorized 
amplitudes, are constrained by the heavy quark effective theory (HQET) 0. Our model 
explicitly satisfies all these constraints and allows the determination of the corrections in 
the inverse powers of the heavy quark masses up to the second order 0. In the quoted 
paper we have determined the Isgur-Wise function and the subleading form factors in the 
hole kinematical range accessible in B —> D transitions. We shall use these functions 
here to evolve B —> D transition form factors from the point of zero recoil of the final D 
meson to the values of q 2 ps rrij , where m/ is a mass of the final light meson. The form 
factors of the heavy-to-light transitions have been calculated in our model at the point 
of maximum recoil of the final light meson using the expansion in inverse powers of the 
heavy b quark mass from the initial B meson and in inverse powers of the large (~ mj,/ 2) 
recoil momentum of the final light meson . We have also determined the q 2 dependence 
of the form factors near this kinematical point. Thus we can calculate the heavy-to- 
light form factors which are necessary for the determination of energetic nonleptonic 
decay amplitudes. This combination of the methods of heavy quark expansion and the 
relativistic quark model increases the reliability of our predictions. The comparison of 
the results with the available experimental data will be the test of factorization. 

The paper is organized as follows. In Sec. 2 we present the expressions for nonleptonic 
decay amplitudes in the factorization approximation. The relativistic quark model is 
described in Sec. 3. In Sec. 4 the heavy-to-heavy transition form factors are discussed. 
The heavy-to-light transition form factors are presented in Sec. 5. Section 6 contain our 
results for the branching ratios of energetic nonleptonic B decays and their discussion. 
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Our conclusions are given in Sec. 7. 


2 NONLEPTONIC DECAY AMPLITUDES AND FAC¬ 
TORIZATION 


In the standard model B decays are described by the effective Hamiltonian, obtained by 
integrating out the heavy W-boson and top quark and applying the operator product 
expansion. For the case of b —> c, u transitions, 


^cb 


He,, = % V ci [cMOf + cMO. 

+% Vel[cMOf + cMOf 


+ 


( 1 ) 


where 14,- are the corresponding Cabibbo-Kobayashi-Maskawa (CKM) matrix elements. 
The Wilson coefficients ci )2 (/f are evaluated perturbatively at the W scale and then 
they are evolved down to the renormalization scale /i fa mb by the renormalization- 
group equations. The ellipsis denotes the penguin operators, which Wilson coefficients 

0 


are numerically much smaller then Ci )2 
are given by 


The local four-quark operators Oj and 0 2 


Of = [(du) V -A + {sc) V - A \(qb)v-A, 

Of = {qu) V - A {db) V - A + {qc) V - A {sb) V -A, q=(u,c), (2) 


where the rotated antiquark fields are 

d Vud,d ~b s V cd d -b 


(3) 


and for the hadronic current the following notation is used 

(qq')v-A = 1 - 75 )q = Jtx- 

The factorization approach to two-body nonleptonic decays B —» Mi(qiq' 1 )M 2 (q 2 q' 2 ) 
implies that the decay amplitude can be approximated by the product of one-particle 
matrix elements: 


(M\ M 2 1 H ef f \B) = ^^V r ^[a 1 (//)(M 1 |(g 1 6)v-A|S)(M 2 |(g 2 ^)v-A|0^ 

+a 2 (//)(M 2 |(g26)v- J 4|5)(Mi|(gi^ 1 ) v _ A |0)], (4) 

where ^ ^ 

= ci(/i) + — c 2 (/r), a 2 (ii) = c 2 (» + — cf/r), (5) 

N c is the number of colors (N c = 3). 
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In the general case, the renormalization point (/i) dependence of the product of current 
operator matrix elements does not cancel the /j dependence of Oj(p) or Cj(/i) [j3], |4j]. Thus 
non-factorizable contributions to ([|) must be present in order to make the physical am¬ 
plitudes independent from the renormalization scale fi. However, as it is shown in [3|], in 
the case of the production of an energetic light meson or meson resonance it is posible to 
justify the factorization approximation and the right-hand side of ([]]) is scale-independent. 
Thus we limit our analysis of nonleptonic decays to consideration of decays with at least 
one energetic meson in the final state (such as B —> D ly *\{p) and B —> 7r(p)7r(p)). 

Before proceeding further, let us additionally note that in writing eq. (f|) we discarded 
the contribution of the colour-octet currents which emerge after the Fierz transformation 
of colour singlet operators (|2j). Clearly these currents violate factorization since they 
cannot provide transitions to the vacuum state. We also neglected the so-called W- 
exchange and annihilation diagrams. In the limit M w —> oo they are connected by the 
Fierz transformation and are doubly suppressed by the kinematic factor of order (' m 2 D /m 2 B ) 
and then dynamicly by the decreasing form factor F D77 (q 2 = m 2 B ) with F Dn (0) = 1 (see 
ref. [Q for details). 

The coefficients (|3]) have been calculated at p ~ mb in the leading logarithmic approx- 
as well as beyond the leading logarithmic approximation H. The result of 


imation 
ref.fffl is 
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«! = 1.01 ± 0.02 and a 2 = 0.20 ± 0.05, 


( 6 ) 


which is close to the result of fitting experimental data 


ai = 1.03 ±0.04 ±0.06 and a 2 = 0.23 ± 0.01 ± 0.01. (7) 

However the a 2 prediction (j5|) is renormalization scheme dependent [[|. 

The matrix element of the current J between the vacuum and final pseudoscalar (P) 
or vector (H) meson states is parametrized by the decay constants fpy 

{P\q r Yn r Y5(l , \Q} = ifpPp, 

(V|?7m^I°) = e P m vfv- (8) 


The matrix elements of the weak current J between meson states have the covariant 
decomposition |l|]: 


(P(p')\<niJ>\B(p)) = 


{V{p')\qinb\B(p)) 


C P + p) 


m 2 B — m 


m B — m 




Fi{q 2 


r 

2 V(q 2 ) . 


q^F 0 (q 2 )] 


m B ± m v 


*IS T fcr 

P P ; 


(9) 

( 10 ) 


{V(p')\qiplsb\B(p)) = [m B ± m v )e* Ai(q 2 ) - Mil {e*q){p + p)^ 

m B + my 
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( 11 ) 


-2 mv^-p-q^A^q 2 ) + 2m v ^^-q^A 0 (q 2 ), 

q z q z 

where q — p — p' and e is a polarization vector of the vector meson. The form factor 
As(q 2 ) is the linear combination 

A I 2\ m B + m v A , bbb - m v 2 
Mq > = ] ~ >' (12) 

and in order to cancel the poles at q 2 = 0, it is necessary to require 

F 1 (0)=F 0 (0), 4»(0) = A 0 (0). (13) 

We calculate the corresponding form factors in the framework of the relativistic quark 
model. 


3 RELATIVISTIC QUARK MODEL 


In the quasipotential approach a meson is described by the wave function of the bound 


quark-antiquark state, which satisfies the quasipotential equation |T3| of the Schrodinger 
type 0 


® *“ (p) = 


V 2 qiR 2 Hr 

where the relativistic reduced mass is 


! j r — 


M 4 


m 


m 


2\2 


4 M 3 


and 


b 2 (M ) = 


[M 2 - (m a + m h ) 2 } \M 2 - (m a - m b ) 2 ] 
4 M 2 


(14) 


(15) 


(16) 


Here m a ^ are quark masses; M is the meson mass and p is the relative momentum of 
quarks. While constructing the kernel of this equation V (p, q; M) — the quasipotential of 
quark-antiquark interaction — we have assumed that the effective interaction is the sum 
of the one-gluon exchange term with the mixture of long-range vector and scalar linear 
confining potentials. We have also assumed that the vector confining potential contains 
the Pauli interaction. The quasipotential is defined by E3 : 


V (p, q, Af) = uJp)u b (-p){-a s D^( k)j^ 

+KLt(k)r;m + v^(k)}«.(9)« l (-9), 


(17) 


where as is the QCD coupling constant, D fiu is the gluon propagator; and u(p) are 
the Dirac matrices and spinors and k = p — q. The effective long-range vector vertex is 
given by 

U(k) = 7m + (18) 
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where k is the Pauli interaction constant. Vector and scalar confining potentials in the 
nonrelativistic limit reduce to 


KL,(r) 

rJ„,M 


(1 - e)(Ar + B ), 
e(Ar + B ), 


(19) 


reproducing V^ c °“ r f el (r) = V^ nf + V^ mf = Ar + B, where e is the mixing coefficient. The 
explicit expression for the quasipotential with the account of relativistic corrections of 
order v 2 /c 2 can be found in ref.|T5[. All the parameters of our model like quark masses, 
parameters of linear confining potential A and B, mixing coefficient e and anomalous 


chromomagnetic quark moment hi were fixed from the analysis of meson masses [|Pj and 
radiative decays [[16]. The quark masses = 4.88 GeV; m c = 1.55 GeV; m s = 0.50 GeV; 
m Ui d = 0.33 GeV and parameters of the linear potential A = 0.18 GeV 2 and B = —0.30 
GeV have standard values for quark models. The value of the mixing coefficient of vector 
and scalar confining potentials £ = —1 has been chosen from the consideration of the 
heavy quark expansion || and meson radiative decays El , which are very sensitive to the 
Lorentz-structure of the confining potential: the resulting leading relativistic corrections 
coming from vector and scalar potentials have opposite signs for the radiative Ml-decays 
» Finally the universal Pauli interaction constant hi = —1 has been fixed from the 
analysis of the fine splitting of heavy quarkonia 3 Pj- states Jl^ . 

The meson wave functions in the rest frame have been calculated by numerical solution 
of the quasipotential equation (|T|) []I7| . However, it is more convenient to use analytical 
expressions for meson wave functions. The examination of numerical results for the ground 
state wave functions of mesons containing at least one light quark has shown that they 
can be well approximated in the meson rest frame by the Gaussian functions 


*" (p)= (S 7 exp {-$. 


( 20 ) 


M, 


with the deviation less than 5%. 
The parameters are 


/ 3 b = 0.41 GeV; 0 D = 0.38 GeV; (3 Ds = 0.44 GeV; (3 <f>) = 0.31 GeV. 


The matrix element of the local current J between bound states in the quasipotential 
method has the form 
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<MU(0)|M) = 


d 3 p d 3 q 
(2vr) 6 


*j/'(p)U(p-q) 4, M(q). 


( 21 ) 


where M(M') is the initial (final) meson, r M (p,q) is the two-particle vertex function and 
are the meson wave functions projected onto the positive energy states of quarks 
and boosted to the moving reference frame. 
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This relation is valid for the general structure of the current = Q'G^Q, where G M 
can be an arbitrary combination of Dirac matrices. The contributions to T come from 
Figs. 1(a) and 1(b). Note that the contribution T*A is the consequence of the projection 
onto the positive-energy states. The form of the relativistic corrections resulting from the 
vertex function T*A is explicitly dependent on the Lorentz-structure of gg-interaction. 

The general structure of the current matrix element m is rather complicated, because 
it is necessary to integrate both with respect to d 3 p and d 3 q. The h-function in the 
expression for the vertex function T (1 l permits to perform one of these integrations. As 
a result the contribution of T^ to the current matrix element has the usual structure 
and can be calculated without any expansion, if the wave functions of initial and final 
meson are known. The situation with the contribution of T® is different. Here instead 
of the (5-function we have a complicated structure, containing the potential of the qq- 
interaction in a meson. Thus, in general case, we cannot perform one of the integrations 
in the contribution of T^ to the matrix element ©■ Therefore, it is necessary to use 
some additional considerations. The main idea is to expand the vertex function T^ in 
such a way that it will be possible to use the quasipotential equation (|T|) in order to 
perform one of the integrations in the current matrix element Q2Tp. The realization of 
such expansion differs for the cases of heavy-to-heavy (B —> anc [ heavy-to-light 

(B —> 7r(p)) transitions. 

4 b —> DECAY FORM FACTORS 

In the case of the heavy-to-heavy (B —> D^) meson decays we have two natural expansion 
parameters, which are the heavy quark masses (mj and m c ) in the initial and final meson. 
The most convenient point for the expansion of vertex function T (2) in inverse powers of 
the heavy quark masses is the point of zero recoil of the final D meson, where A = 0 
(A = p b — p dm). If is easy to see that D 2 ^ contributes to the current matrix element 
at first order of the 1 /rriQ expansion. We limit our analysis to the consideration of the 
terms up to the second order. After the expansion we perform the integrations in the 
contribution of T® to the decay matrix element. As a result we get the expression for the 
current matrix element, which contains the ordinary mean values between meson wave 
functions and can be easily calculated numerically. The results of such calculation are 
given in comparison with the predictions of HQET [I7j in [|]. Our model satisfies all 
the constraints imposed on the form factors by heavy quark symmetries and allows the 
determination of the Isgur-Wise and subleading form factors j8|. 

The g 2 dependence of form factors at leading order of the 1/rriQ expansion is given by 

m(q 2 ) = R F ° (q2 J = R‘V(q 2 ) 

m B +m D 

= R*A 0 (q 2 ) = R* Mq2 J = R* A 2 (q 2 ) = £(w), (22) 

m B +m D * 
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where 


R (*) = 


2 y/m B m D (*) 


w = 


m B + m D (*) 
m B + m DM - Q 2 


2 m B m D (*) 


The Isgur-Wise function in onr model is 


= 


1\ w — 1 


— exp ^-^2, - 2Jw+1 


( 23 ) 


with the slope parameter p 2 ~ 1.02, which is in accordance with the recent CLEO II 
measurement [[E| p 2 = 1.01 ± 0.15 ± 0.09. 

At the Erst order of the 1/rriQ expansion four additional independent form factors 
arise in HQET We have determined these subleading form factors in the framework 
of our model HI 


60) 
XiM 
X 2 (w) 
XsM 


= (A - m q ) (1 + 


2w- 1 
3 w + 1 




1 A 

32 w + 1 
1 - w — 1 


= —A 


16 w + 1 


fM, 


(24) 


where the HQET parameter A — M — rriQ in our model is equal to the mean value of the 
light quark energy in the heavy meson A = (e q ) ~ 0.54 GeV. 

We have also calculated the second order power corrections at the point of zero recoil of 
the final meson §|. The obtained structure of the second order corrections is in accord with 

IT*) transition form 


predictions of HQET [pTS 
factors at q 2 = q 2 


As a result we got the values of the B 
up to the second order terms. The higher order terms of the l/m,Q 
expansion are negligibly small. However, for the consideration of the energetic nonleptonic 
decays of B mesons we need the form factor values at q 2 = m 2 ~ 0 (/ = tt, p... is a final 
light particle). Thus it is necessary to evolve the form factors from g 2 iax to q 2 ~ 0. The 
corresponding w range is not very wide (from 1 to ~ 1.6). For such w values the form 
factors are dominated by the Isgur-Wise function |20|. Some small contributions may 
arise from subleading form factors. The higher order terms give very small corrections 
2|. Therefore we combine the universal Isgur-Wise q 2 dependence of form factors (|22|) 
with the subleading symmetry breaking corrections (|24|). The resulting form factor w 
dependence is shown in Figs. 2,3. 

The values of form factors at q 2 = 0 are 


iA(0) = F 0 (0) = 0.63, F(0) = 0.79, 

A 0 (0) = 0.63, A^O) = 0.62, A 2 (0) = 0.61. 


( 25 ) 





















5 B —* 7r(p) DECAY FORM FACTORS 


In the case of heavy-to-light decays the final meson contains only light quarks (u, d ). 
Thus, in contrast to the heavy-to-heavy transitions, we cannot expand matrix elements 
in inverse powers of the final quark mass. The expansion of T^ only in inverse powers 
of the initial heavy quark mass at A = 0 does not solve the problem. However, the final 
light meson has the large recoil momentum, in comparison with its mass, almost in the 
whole kinematical range. At the point of maximum recoil of the final light meson the 
large value of recoil momentum |A max | ~ m\,j 2 allows for the expansion of decay matrix 
element in 1 /mj. The contributions to this expansion come both from the inverse powers 
of heavy mb from the initial B meson and from inverse powers of the recoil momentum 
| A max | of the final light 7r(p) meson. In ref. [H we carried out this expansion up to the 


second order and performed one of the integrations in the current matrix element (|2l|) 
using the quasipotential equation as in the case of a heavy final meson. As a result we 
again get the expression for the current matrix element, which contains only the ordinary 
mean values between meson wave functions, but in this case at the point of maximum 
recoil of the final light meson. 

The found values of B —> tt( p) form factors at the point of maximum recoil (g 2 = 0) 
are 


Fi(0) = F 0 ( 0) = 0.21, V(0) = 0.29, 


4,(0) 

= 0.18, A(0) = 0.27, A 2 (0) = 0.30. 

(26) 

The g 2 behaviour of the heavy-to-light form factors near A max 

(corresponding to 

g 2 = 0) is given by 0 



Fi{q 2 ) 

M b + Mn ~ f A , A 2 ^ 

2 ( )^ + ( max ’ 

(27) 

Fo(q 2 ) 

M B + M p 2 1 max), 

(28) 

Mq 2 ) 

2 X IM b M p i 

= m b + m~A + w ^ m ^’ 

(29) 

A ),2 (q 2 ) 

= 2 ( a^j, 

2 y M B M p 

(30) 

V{q 2 ) 

Mb + M p ~ u 2 ^ 

= „ nrj , f £M v (A max ), 

2 J Mb Mp 

(31) 


where we have introduced the function || 

I 2 [ A 2 w — l\ 

= (32) 
9 ft 2 

V = , A » 0.53 GeV. 

Ab + Ptx (p) 
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Equation (^) reduces to the Isgur-Wise function (^) in the limit of infinitely heavy 
quarks in the initial and final mesons. 

It is important to note that the form factors Ai and F 0 in (^), (^8|) have a different 
q 2 dependence than the other form factors (pi7|), ([¥]), ([IT]). In quark models one usually 
assumes a pole |21|] or exponential p2[ q 2 behaviour for all form factors. However, the 
recent QCD sum rule analysis indicates that the form factor A\ has a q 2 dependence 


24, 25 


which is different from other form factors | 

The extrapolation of the q 2 dependence (|28|) -(|3T1) to all values of q 2 (or w), accessible 
in B —> 7r(p) transitions, introduces rather large uncertainties, because w varies in a broad 
kinematical range (from 1 to ~ 19 in B —* tv and from 1 to ~ 3.5 in B —> p). However, 
the w values for B —> tv (p) form factors which are really necessary for the consideration 
of energetic nonleptonic decays are limited to a rather small interval near w max ( q 2 = 0). 
Thus, the application of formulas (|28|)-(|3T|) in this region is rather reliable. We show the 
w dependence of B — > tv( p) decay form factors in Figs. 4,5. 


6 RESULTS AND DISCUSSION 


In the factorization approximation one can distinguish three classes of B meson nonlep¬ 
tonic decays (see Fig. 6) [[[]: the ‘class I’ transitions, such as B° —> , where only 

the term with cq in (|4|) contributes (i.e. both mesons are produced by charged currents); 
‘class II’ transitions, such as B° —> M^M®, where only the term with a 2 in ([!]) contributes 
(i.e. both mesons are produced by neutral currents) and ‘class III' transitions, such as 
B~ —> , where both terms can contribute coherently. 

The results of the calculation of the nonleptonic branching ratios for B —> and 

B —> 7r(p) transitions are given in Tables 1 and 2 in comparison with other model pre¬ 
dictions 0-0 and experimental data. The B —> /}(*)£)(*) decay branching ratios are 
presented for completeness. 

We see that our results for the ‘class I" nonleptonic B —> D^Tv{p) decays are close 
to the improved BSW model predictions f26|, while our results for the ‘class II' and a 2 
contributions to ‘class III’ decays are smaller than those of [^]. These contributions come 
from B —> Tv(p) transition form factors, which have a different q 2 behaviour in our and the 
BSW models. The BSW model assumes universal q 2 dependence of all B —> iv(p) form 
factors. As already mentioned in our model we find the A\ and Fq form factors to have a 
different q 2 dependence than that of the other form factors (see (|2E|)-(|3T|) and Figs. 4,5). 
The form factor F 0 in our model decreases with the growing of q 2 (decreasing of w) in 
the kinematical range of interest for energetic nonleptonic decays (see Fig. 4). Note that 
our value for B —> tv form factors at q 2 = 0 is approximately 1.5 times less than that of 
BSW, while the values for B —> p form factors are close in both models. 

Our predictions for the branching ratios of B —* D^M nonleptonic decays presented 
in Table 1 agree with experimental data within errors. Thus we can conclude that factor¬ 
ization works rather well for ‘class I’ and ‘class III’ decays B —> D^tv(p). However, an 
improvement of experimental accuracy is needed to make a definite conclusion. It will be 
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Table 1: Predicted branching ratios for B —> D^M nonlcptonic decays in terms of cq 
and a 2 (in %). Onr model branching ratios are quoted for values of eq = 1.05 and 
02 = 0.25 in comparison with experimental data (in %). We use the values \V c b\ = 0.038 
and fo = /d* = 0.220 GeV, fo s = fD* = 0.260 GeV, f ai = 0.205 GeV [12] for our 


estimates. 


Decay 

our result 

26] 

our result 

experiment [29] 

5 U -*■ 

D + n 

0.29aj 

0.264a? 

0.32 

0.31 ± 0.04 ± 0.02 

B° -> 

D+p~ 

0.79af 

0.621a? 

0.87 

0.84 ±0.16 ±0.05 

B° 

D*+7T~ 

0.26af 

0.254a? 

0.28 

0.28 ±0.04 ±0.01 

B° 

D*+p~ 

0.81a? 

0.702a? 

0.88 

0.73 ±0.15 ±0.03 

B° -*■ 

D + a{ 

0.78a? 

0.673a? 

0.86 

0.60 ±0.22 ±0.24 

B° -*■ 

D*+a{ 

0.92a? 

0.970a? 

1.02 

1.27 ±0.30 ±0.05 

B° 

D + D s 

1.37a? 

1.213a? 

1.51 

0.74 ±0.22 ±0.18 

B° 

D+D*~ 

0.685a? 

0.859a? 

0.75 

1.14 ±0.42 ±0.28 

B° 

D*+D~ 

0.82a? 

0.824a? 

0.90 

0.94 ±0.24 ±0.23 

B° -»• 

D*+D*~ 

2.50a? 

2.203a? 

2.75 

2.00 ±0.54 ±0.49 

B° -»• 

D°ir° 

0.058a? 

0.20a? 

0.0036 

< 0.048 

B° 

D*°ir° 

0.056a? 

0.21a? 

0.0035 

< 0.097 

B° 

D°p° 

0.053a? 

0.14a? 

0.0033 

< 0.055 

B° -»• 

D*°p° 

0.156a? 

0.22a? 

0.0098 

< 0.117 

-* 

D°ir~ 

0.29(ai + 0.64a 2 ) 2 

0.265(ai + 1.230a 2 ) 2 

0.43 

0.50 ± 0.05 ± 0.02 

B~ -» 

D*° 7T- 

0.27(ai + 0.69a 2 ) 2 

0.255(ai + 1.292a 2 ) 2 

0.40 

0.52 ± 0.08 ± 0.02 

-» 

D°p~ 

0.81(ai +0.36a 2 ) 2 

0.622(ai +0.662a 2 ) 2 

1.06 

1.37 ±0.18 ±0.05 

B~ -> 

D*°p~ 

0.83(a? + 0.39a? 

+ 1.15aia 2 ) 

0.703(a? + 0.635a? 

+ 1.487aia 2 ) 

1.17 

1.51 ±0.30 ±0.06 

B~ -* 

D°D~ 

1.40a? 

1.215a? 

1.55 

1.36 ±0.28 ±0.33 

B~ -* 

D°D*~ 

0.70a? 

0.862a? 

0.77 

0.94 ±0.31 ± 0.23 

B~ -» 

D*°D~ 

0.84a? 

0.828a? 

0.92 

1.18 ±0.36 ±0.29 

B~ -» 

D*°D*- 

2.56a? 

2.206a? 

2.80 

2.70 ±0.81 ±0.66 


very interesting to measure the ‘class If’ decay B° —> D^°tt(p) 0 branching ratios. Such 
measurement will be the test of factorization for ‘class 11’ nonleptonic decays and will 
help to constrain the w (or q 2 ) dependence of B —> 7r(p) form factors. 

We also present in Table 1 the predictions for B —> nonleptonic decays, where 

only heavy mesons are present in the final state. The factorization is less justified for such 
decays. However, as we see from Table 1 our predictions, based on the factorization, are 
consistent with the experimental data for these decays too. 

For the branching ratios of B —> 7i(p)M nonleptonic decays presented in Table 2 there 
are only experimental upper limits at present. The measurement of these decays will allow 
the determination of the CKM matrix element \V u b\, which is poorly studied. The closest 
experimental upper limit to the theoretical predictions is for the decay B° 7r + 7i~. ft is 
approximately two times larger than our prediction and is very close to the result of El- 
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Table 2: Predicted branching ratios for B —> 7r(p)M nonlcptonic decays. We use the 
experimental values for / f p [12] and the value of |Kb| = 0.0052 [9] for our model 
estimates. All numbers are branching ratios xlO 5 . 


Decay 

B° —■> 7r + 7T _ 

B° -> vr+p' 
B° —■> p + TT~ 
B° -> 7T ± p^ 
5° -> p+p- 

5° -> 7T 0 7r 0 

B° —■> 7r°p° 
5° -> p°p° 

5° -> D~n + 
5 ° -> D*~7T^~ 
B° D~p + 
B° -»■ D*-p+ 

5” -> 7T°7r- 
—■> 7T°p — 
5” —■> p°7T~ 
-> p°p" 

B~ D S -7T° 

D*-7T° 

5” D s -p° 

b- - or/ 


our result 


0.331 |K* 

2 a 2 

0.857| Kb 

| 2 a 2 

0.234 Kb 

| 2 a 2 

1.09|Kb| 

2 a 2 

0.794 Kb 

| 2 a 2 

0.17 Kb 

2 a 2 

0.54 Kb 

2 a 2 

0.39|1/1 

to to 

0.285 Kb 

| 2 a 2 

i.06|K & | 

2 a 2 

0.269 Kb 

| 2 a 2 

2.25 Kb 

2 a 2 


0.169|Kb| 2 (ai + CI 2) 2 
0.438|Kfe| 2 (ai + 0.52a 2 ) 2 
0.120|Kfe| 2 (ai + 1.95a 2 ) 2 
0.411|K b | 2 (ai + a 2 ) 2 
0.146|Kb| 2 a? 
0.545|Kfe| 2 a? 

0.138 Kb 2 a? 

1 . 16 |K b | 2 a? 


our result []27| [28|] experiment UL 


0.99 

1.8 


< 2.0 

2.55 

4.8 



0.70 

0.4 



3.25 

5.2 


< 8.8 

2.36 

1.3 


< 220 

0.028 

0.06 


< 0.91 

0.092 

0.14 


< 2.4 

0.067 

0.05 


< 28 

0.85 

8.1 

1.9 

< 28 

3.1 

6.1 

2.7 

< 50 

0.80 

1.2 

1.0 

< 70 

6.7 

4.5 

5.4 

< 80 

0.78 

1.4 


< 1.7 

1.7 

2.7 


< 7.7 

0.77 

0.7 


< 4.3 

1.8 

1.1 


< 100 

0.44 

3.9 

1.8 

< 20 

1.6 

3.0 

1.3 

< 33 

0.41 

0.6 

0.5 

< 40 

3.4 

2.4 

2.8 

< 50 



From this upper limit on B(B° n + i r ) we get the limit on | Kb I in our model 

|Kb| < 7.4 x 10 -3 , 

which is close to the value previously found from semileptonic B —> tt(p)Iu decays 0: 

\V uh \ = (5.2 ±1.3 ±0.5) x 10~ 3 . 

7 CONCLUSIONS 

In this paper we have calculated the branching ratios of the energetic exclusive nonlcptonic 
decays of B mesons on the basis of the factorization approximation. In particular the form 
factors of B —> and B —> 7r(p) transitions have been evaluated using the relativistic 

quark model and the heavy quark expansion. Such expansion has been carried out up 
to the second order in the heavy quark masses. Finally the momentum dependence of 
leading and subleading terms of this expansion has been used for the determination of 
the heavy-to-heavy transition form factors at q 2 = m 2 , where / is a final light meson. 


12 















The overall agreement of our predictions for two-body nonleptonic decays of B mesons 
with the existing experimental data ||29|. |3(J shows that the factorization approximation 
works sufficiently well in the framework of our model. From another side, if the fac¬ 
torization hypothesis is taken for granted, the aforementioned agreement confirms the 
selfconsistency of our approach, which incorporates our previously obtained results for 
semileptonic and leptonic decays of heavy mesons. In particular it would be quite in¬ 
tersting to test the specific q 2 behaviour of the heavy-to-light transition form factors Fq 
and A\ predicted by our model. Another important problem is the determination of the 
coefficients a± and a 2 via c\ and c 2 directly from QCD. As it has been discussed already 
in Sec. 2 their values found in ref. || are close to those obtained from fitting experi¬ 
mental data, though one should mention that c 2 is rather unstable with respect to the 
renormalization scheme and scale. Nevertheless that means from our point of view that 
the factorization hypothesis for energetic B meson decays has more or less firm grounds 
within QCD, at least for the ‘class I’ decays. 

The situation is essentially different for D meson nonleptonic decays. In this case 
the best fit to the experimental data yields oq = 1.26 ± 0.10 and a 2 = —0.51 ± 0.10 
0.0. Meanwhile QCD predicts [0] (the instability of the coefficient here is even stronger) 


Ci = 1.31 ± 0.19 and c 2 = —0.50 ± 0.30, which can only be consistent with the result of 
fitting data if one drops the 1 /N c terms in eq. (0) and puts cq ~ Ci and a 2 ~ c 2 . Clearly, 
such an assumption would give a completely wrong result for B decays, namely a negative 
sign of a 2 , which is ruled out by the experimental data. This result could indicate that 
the factorization approach in D decays is insufficient and non-factorizable contributions 
there are large. In other words the D meson is possibly not heavy enough compared to 
the B meson. The results of refs. (J3J seem to point in the same direction. 

The successful description of nonleptonic two-body decays of B mesons makes the 
present approach appealing for the further consideration of B s c meson decays. Work in 
this direction is in progress. 
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FIGURE CAPTIONS 


FIG. 1. (a) The lowest order vertex function (b) The vertex function T® with the 
account of the quark interaction. The dashed line corresponds to the effective potential 
(IT). The bold line denotes the negative-energy part of the quark propagator. 

FIG. 2. The w dependence of form factors of B —» D transitions. 

FIG. 3. The w dependence of form factors of B —> D* transitions. 

FIG. 4. The w dependence of form factors of B —> n transitions in the kinematical region 
of interest for the energetic nonleptonic decays. 

FIG. 5. The w dependence of form factors of B —> p transitions in the kinematical 
of interest for the energetic nonleptonic decays. 

FIG. 6. Quark diagrams for two-body nonleptonic B decays: (a) ‘class I" decays 
Mi M 2 ; (b) ‘class II’ decays B° —> (c) ‘class III’ decays B~ —> AT) 1 ATT. 


region 
B° -> 
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